Anisotropic in-plane resistivity in the nematic phase of the iron pnictides 
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We show that the interference between scattering by impurities and by critical spin fluctuations 
gives rise to anisotropic transport in the Ising-nematic state of the iron pnictides. The effect is closely 
related to the non- Fermi liquid behavior of the resistivity near an antiferromagnetic quantum critical 
point. Our theory not only explains the observed sign of the resistivity anisotropy Ap in electron 
doped systems, but also predicts a sign change of Ap upon sufficient hole doping. Furthermore, our 
model naturally addresses the changes in Ap upon sample annealing and alkaline-earth substitution. 



In many materials, anisotropic properties are related 
to the underlying crystalline structure. However, when 
correlations are present, the electronic states can them- 
selves become anisotropic fH . Recent transport measure- 
ments in detwinned crystals of the iron pnictide com- 
pomids ^i?Fe2As2 {AE = Ba, Ca, Sr) found an in-plane 
anisotropy that cannot be attributed only to lattice dis- 
tortions, unveiling an anisotropic electronic state 
Its existence is also supported by the observations of an 
orbital polarization O — rij-z—Uyz ofd^z and dyz Fe states 
in angle-resolved photoemission spectroscopy (ARPES) 
4 1, local anisotropics in scanning tunneling microscopy 
5|, and anisotropics in the optical spectrum 
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One candidate for such unconventional electronic state 
is the Ising-nematic order, which emerges from the com- 
bination of magnetic fluctuations and frustration 
In a strong-coupling approach, frustration is promoted 
by the competing Ji — J2 exchange interactions, while 
in an weak-coupling approach, it follows from the degen- 
eracy of the magnetic ground state due to the nesting 
properties of the Fermi surface (FS). Indeed, the iron 
pnictides support two magnetic instabilities with order 
parameters Ai, A2 corresponding to the in-plane order- 
ing vectors Qi — (tt, 0) , Q2 — (0, tt), respectively. The 
electronic structure [lj,ll5|, as well as the coupling to the 



lattice jl3|, give rise to the coupling ^Af in the free 



energy, with g > 0. As a result, a discrete Ising-nematic 
degree of freedom tp oc (Af ) — (A|) emerges, which labels 
the two degenerate states Ai 0, A2 = and Ai = 0, 
A2 7^ 0. Most importantly, ip is able to order at a temper- 
ature Tg above the onset of antiferromagnetism (AFM) 
at Tn jl6l|, breaking the tetragonal symmetry already in 
the paramagnetic (PM) phase - hence the term nematic. 

By symmetry, Ising-nematic order induces the orbital 
polarization O and a shear distortion Eg , driving a struc- 
tural transition from the tetragonal (Tet) to the or- 
thorhombic (Ort) phase at Tg. It explains why Tg and 
Tjv track each other closely across every pnictide phase 
diagram, even inside the superconducting (SC) dome 
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Recently, it has been shown that nematic fluc- 
tuations can explain the dramatic softening of the lattice 
observed experimentally in the Tet phase |13| . It remains 



Figure 1: (a) Schematic temperature-doping {T,x) phase di- 
agram in the absence of SC. The yellow region corresponds 
to enhanced nematic fluctuations, as seen by elastic modu- 
lus measurements [Tsl. (b) Temperature dependence of the 
correlation lengths ^ associated with the (tt, 0) and (0, tt) mag- 
netic instabilities above the magnetic QCP. 



an important issue whether this scenario is also able 
to address the observed in-plane resistivity anisotropy 
Ap = pb — Pa, where b (a) refers to the direction with 
shorter (longer) lattice constant of the Ort state. The in- 
duced orbital polarization O has been alternatively pro- 
posed to explain such anisotropy, but different authors 
disagree on the suitability of this proposal (l9l . [20| . 

In this paper we demonstrate that the observed 
anisotropic transport above T^v can be understood within 
the Ising-nematic picture. Our model explains in a uni- 
fied way several experimental observations, such as the 
sign of the anisotropy 0, and the changes in Ap upon 
sample annealing [q| or alkaline-earth substitution |2l] |. 
We also predict a sign change of Ap for hole-doped sys- 
tems, thus addressing the peculiar observations in K- 
doped compounds [22, 23]. Our results support the view 
that nematic order is the driving force for the structural 
transition and the orbital polarization in the pnictides. 
Note that our theory does not address the anisotropy ob- 
served below T/v, in the AFM state. However, in all ob- 
servations, the sign of Ap is already set above Tn, where 
the slope d{Ap) /dT is the largest [3,|1|. Concentrating 
on the regime above allows for a more transparent 
and less model-dependent understanding of the resistiv- 
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ity anisotropy. 

Below Ts, the onset of a finite Ising-nematic order pa- 
rameter 7^ gives rise to two sources of anisotropy: in 
the electronic structure itself and in the spectrum of mag- 
netic fluctuations. The former is rather weak, due to the 
smallness of the shear distortion Eg — {a — b) / {a + b) 
- for instance, in optimally doped Ba (Fei_xCox)2 As2, 
sa 5 X 10~^ 0, Il3|- On the other hand, the latter 
is strong, since only fluctuations associated with one of 
the magnetic instabilities diverge. The impact of this 
anisotropy on the transport depends on the contribu- 
tion of the scattering of electrons by spin fluctuations 
around the hot spots of the Fermi surface, defined via 
ex.k — £A'.k+Qi, with ex,k denoting the dispersion of 
band A. In ultra-clean systems, this scattering channel is 
short-circuited by the contribution from the other regions 
of the FS [23|. However, in the presence of significant 
impurity scattering, the scattering by spin fluctuations 
becomes important, leading to non- Fermi liquid behav- 
ior [2^ and, in addition, to anisotropic transport. To 
see this, we follow Rosch [2^ and solve the Boltzmann 
equation (BE) for scattering by impurities and spin fluc- 
tuations, which yields features not captured by the re- 
laxation time approximation [2^. In the regime where 
impurity scattering dominates, we obtain the resistivity 
along the a {a = x, y, z) direction: 

kk'.AA' 

With V^^, ^ /o (1 - fl^ t^,%,. Here, k is the elec- 

tron momentum, A is the band index, ^ is the Fermi- 
Dirac distribution function, and pi,np is the (isotropic) 

residual resistivity, i^f^^.k' ^ 9sUin il,^) 

is the spin-fluctuations collision term, where uj — e\^i^ — 
and q = k — k' are, respectively, the transferred 
energy and momentum, and Xi (Qi '-^) is the dynamic sus- 
ceptibility associated with the ordering vector Q^. g^f\ 
is the scattering amplitude and ^ is the deviation of 
the electron distribution function from equilibrium due 
to scattering by impurities only, as a consequence of an 
electric field applied along the a direction. It only de- 
pends on the impurity scattering potential V^^p and on 
the shape of the FS (see supplementary material for more 
details). 

The resistivity anisotropy data in Ba (Fei_xCox)2 As2 
show that the temperature-dependent contribution is 
smaller than thejesidual resistivity, particularly close to 
optimal doping [2]. Therefore, we can calculate paa us- 
ing Eq. whose inputs are the band dispersions £A.k, 
which determine the impurity-only solution ^^'^ 
the dynamic susceptibilities Xi(q, w), which determine 
the spin- fluctuation scattering amplitude t^^^^ j^/ . In the 
Ising-nematic phase, the latter are: 



ky(*)» 




Figure 2: Electrons at the active and passive hot spots (full 
and empty circles, respectively) are scattered with different 
amplitudes by spin fluctuations associated with the (tt, 0) and 
(0, tt) ordering vectors. The projection of the hot spots Fermi 
velocities (arrows) along the (a, b) axes determine the sign of 
the low-temperature resistivity anisotropy. 



Xi (q, = , , r> N (2) 

(q + Qi) T <y5 - «^^/7 

where the upper (lower) signs refers to Qi (Q2), 7 
is the Landau damping parameter, and (q + Q^) ~ 
aU-' + olll (1 ± ^) + ^Wv (1 T ry) -I- r;^ cos^ (i^) . This 
phenomenological form follows naturally from the elec- 
tronic structure of the iron pnictides [27j . and is con- 
firmed by inelastic neutron scattering experiments (28| . 
Here, we introduced the in-plane (out-of-plane) tetrag- 
onal lattice constant oo (cq), the magnetic correlation 
length ^, the in-plane momentum anisotropy -q, and the 
out-of-plane anisotropy rjz . The temperature dependence 
of ^ and Lf are calculated within a self-consistent mean- 
field approach [ill. [l2|. yielding the equations: 

- = [Xl (q,W«) - X2 (q,W„)] /istrain, (3) 

q I ^ — ' 

where a;„ are bosonic Matsubara frequencies and 
^strain P, with P denoting the external strain applied 
to detwin the samples. For /istrain = 0, Eqs. [3] have 
a non-trivial solution with (/? 7^ below a temperature 
Ts, describing the PM-Ort (Ising-nematic) phase. This 
model is similar to the one used in Ref. [ISl, which suc- 
cessfully describes the lattice softening above Tg. Since 
^j^^ = ^ if, when > (i.e. a || x and 6 || y) the 
correlation length f(7r,o) associated with the Qi ordering 
vector lengthens while C(o,7r)j associated with Q2, short- 
ens (see Fig. 1). In the figure, we present results for an 
AFM quantum critical point (QCP), C()r^o) = 0) = 0' 
where fluctuations are strongest. 



3 



Substituting the susceptibilities given by ^ in the BE 
solution ([l}, we find that the low-T transport is domi- 
nated by the subset of hot spots connected by the soft 
ordering vector (Qi for (p > 0), e^.k^, = Eamc^+Qi , while 
contributions due to Q2 = (0, tt) are negligible. It follows 
that 

Paa - Amp - - Km.X')' (4) 

where k > is a constant. The exponent C deter- 
mines the T-dependence of the resistivity and is given 
by C = 3/2 (C = 1) for T <Ar]l (T > K-ql), depending 
on whether the spin fluctuation spectrum is two- or three 
dimensional (A is the energy cutoff). While the experi- 
ments of Ref. [2^ on BaFe2 (Asi_xPx)2 materials support 
C = 1, Eq. (HI) is always dominant when compared to the 
Fermi- liquid contribution p — Pimp oc T^. The distinction 
into active and passive hot spots is not restricted to the 
AFM QCP and persists everywhere along the finite tem- 
perature phase boundary where soft and hard magnetic 
fluctuations exist. Thus, the Ising-nematic order causes 
an anisotropy in the spectrum of magnetic fluctuations, 
which induces an anisotropy in the resistivity via scat- 
tering by spin fluctuations (see Fig. 2). To determine the 
amplitude and sign of Ap, all that is left is to calculate 
i.e. the solution of the BE with impurities only, 
which depends solely on the band dispersions e^.k- 

To gain analytical insight, we introduce a simplified 
model with a circular hole pocket at the center of the Bril- 
louin zone and elliptical electron pockets displaced from 
the center by Qi and Q2 ^i,k = £0 — fc^/2TO — p, 

£2.k+Qi = -£o + kl/2mx -I- kl/2my - /i, and £3,k+Q2 = 
—Co -l- k'^/2my + ky/2iTix — /i. Here, Eq is the shift in 
the band energies, rrii are the band masses, and fj,, the 
chemical potential. The calculation of $g ^ for this case 
is straightforward and is presented in the supplementary 
material. In particular, we find that the sign and ampli- 
tude of the anisotropy depends on the projections along 
kx (parallel to a) and ky (parallel to b) of the Fermi ve- 
locities at the hot spots ei.k^^, — £2,kfe,-i-Qi (see Fig. 2). 
When the hot spots are close to a certain axis, there is 
stronger scattering by spin fluctuations when the elec- 
trons move parallel to this axis, and the resistivity is 
larger along this direction. Interestingly, ARPES mea- 
surements in the optimally doped Ba (Fei_xCox)2 As2, 
which is close to a possible AFM QCP covered by the 
SC dome [30], reveal the presence of hot spots close to 
the ky (b) axis [sH . According to our model, this implies 
Pb > Pa, in agreement with transport measurements 

To go beyond the low-T limit in Eq. we numeri- 
cally calculated Ap from Eq. ^ for fixed band structure 
parameters, mx = 1.25m, ruy — 0.83m, and p = O.OSeo. 
For the dynamic susceptibility in Eq. ([2]), based on previ- 
ous works and neutron scattering data [i^, we used 
g/Co w 2 X 10-2, « 0.1, ry w 0.5, 7 w 350 meV, 
and Co « 1 meV, yielding T, = AO K and Tn = 0. We 
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Figure 3: Temperature dependence, above the magnetic QCP, 
of the resistivities along a [pa, green curve) and h {pb, red 
curve), in the absence (dashed lines) and presence (full lines) 
of external strain. The inset shows the resistivity anisotropy 

Ap = Pb - Pa- 

also considered g^f\ sa 2V^^^ and expressed the results 
in terms of the residual resistivity pimp. 

The behavior of pa and pt as function of temperature 
is shown in Fig. 3, where we considered (p > 0, i.e. x 
and y parallel to a and 6, respectively. In accordance 
to our discussion of the low-T limit, we find ph > pa 
for any temperature below T^, a behavior that lingers 
even when one moves away from the QCP. As expected, 
we also find that Ap ~ ip close to Tg [32]. Note that 
the sign of Ap does not depend on the sign of the in- 
plane momentum anisotropy 77 of x in Eq. As 
T — >■ 0, spin fluctuations are suppressed and Ap — >■ 0. In 
the actual data of Ref. [2] , these curves are cut off by the 
onset of SC, which is not considered in our model. 

To make the comparison with experiments more realis- 
tic, we also included the effects of a small strain P = 0.1 
MPa, which is applied to detwin the sample. This strain 
breaks the tetragonal symmetry at any temperature, and 
there is not a well-defined structural transition, since 
(p ~ /ig/rain finite at T*™", where S is the correspond- 
ing Ising critical exponent. Our calculations show that 
the AFM transition is practically unchanged for small 
strain, in agreement with the observations in Ref. 
If the structural transition in the twin sample is strongly 
first order such that ft-g/rlm small compared to the jump 
of If at T*"'", the tail in Ap above T^^'"" becomes negli- 
gible. This is the case in the parent compounds with Ba 
replaced by Ca or Sr, as observed by Ref. [211] . 

Our model also explains recent experiments that found 
a suppression of Ap after annealing the sample fs]. As 
we mentioned earlier, in the ultra-clean limit the contri- 
bution of the hot spots to the resistivity is short-circuited 
by the other regions of the FS [24]. Since the anisotropy 
is governed by scattering processes near the hot spots, 
Ap becomes smaller for cleaner samples. Our results can 
also be rationalized in terms of a T-dependent anisotropic 
dressing of the impurity scattering by spin fluctuations. 



4 




-0.2 -0.1 

hole doping - 



electron doping 



Figure 4: Sign of the low-temperature resistivity anisotropy 



Ap — pb — Pa as function of the ellipticity 



of the elec- 



tron pockets and the chemical potential fi. Note the asym- 
metry between electron and hole doping. 



Via this many-body mechanism, an s-wave scattering 
center acquires effectively an anisotropic cross-section, 
breaking the C4 symmetry. 

Finally, we discuss a general prediction of our model. 
When doping is introduced and the chemical potential 
shifts, the positions of the hot spots change accordingly. 
Without relying on details of Sa.Ic, one expects that by 
hole-doping the samples (i.e. decreasing fi towards nega- 
tive values), the area of the hole pocket increases, whereas 
the area of the electron pocket decreases. Consequently, 
active hot spots, which initially are close to the ky (b) 
axis, will move towards the fcj. (a) axis, and Ap will even- 
tually change its sign. This is illustrated in Fig. 4, where 
we present a phase diagram for the sign of Ap for different 
values of the chemical potential fi and the ellipticity of 
the electron pockets cx mx—niy, keeping the other param- 
eters unchanged. It is interesting to note that measure- 
ments of Ap in hole-doped samples (Bai_xKx) Fe2As2 
found a vanishingly small resistivity anisotropy [22!], as 
expected near the region of the phase diagram of Fig. 4 
where Ap changes sign. Most interestingly, very recent 
data on samples with higher doping revealed a negative 
Ap, in accordance to our prediction (23| . 

In our approach, the resistivity anisotropy above T/v 
is a consequence of anisotropic scattering rates along a 
and b directions. Thus, optical conductivity at low fre- 
quencies is the ideal tool to verify our results. Unfortu- 
nately, as stated in Ref. 01, the experimental conditions 
in this temperature regime need to be refined before def- 
inite conclusions can be drawn. Yet, it is encouraging 
that the available data above T/v on electron-doped sam- 
ples indicate a larger scattering rate along the b direction 

0|, in agreement with our results. We stress that our 
model is suitable for the paramagnetic phase, where the 
sign of the resistivity anisotropy is determined, accord- 
ing to the experimental results on various compounds 
0, Sj HH, [15. Below Tjsf, the spin fluctuations associ- 



ated with the magnetically ordered phase are naturally 
anisotropic, and also give rise to anisotropic scattering. 
Both this anisotropic scattering and the details of the 
anisotropic reconstruction of the FS (20l . 33, 34| will gov- 



ern the temperature evolution of Ap below T/v. 

In summary, we showed that the onset of Ising-nematic 
order leads to an anisotropy in the magnetic fluctuations 
associated with the two magnetic ground states of the 
iron pnictides. As a result, the same physics responsible 
for non-Fermi liquid behavior - the interference of scat- 
tering by spin fluctuations and impurities |2^ - also gives 
rise to anisotropic transport properties. The latter are in 
agreement with several different experimental observa- 
tions, such as the opposite signs of Ap in electron-doped 
and hole-doped materials, as well as the suppression of 
Ap upon sample annealing. 
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Supplemental Material 

The linearized Boltzmann equation determines the electronic distribution function /k,. 



k'.A' 



($ 



4> 



k'.A 



) /k',A' (l /k,A) 



,A.A' 
''k.k' 



(SI) 



Here, <I>k,A is the non-equilibrium distribution function, defined as /^.a = /k a ~ (95^) ^k,A: with ^ denoting 
the Fermi-Dirac distribution at energy Ek.A- k is the momentum, A is the band index, Vk.A = d£\^,\/dk., and n is the 
unit vector parallel to the external electric field. The collision term t^'^, comprises scattering by impurities and spin 
fluctuations [23.l25t: 



^fmp.k.k' = (5?mp'^A,A' + ffimp (1 - '5a,A')) ^ (Sk.A - £k',A') 



imp.k.k 
,A.A' 



sf,k,k' 



2g2 

— J^Cck.A - £k',A')ImX» (k- k',ek,A - £k',A') 



(S2) 



with intraband impurity scattering amplitude gf^.^ 



\=\' 



interband impurity scattering amplitude g^^^ = 



imp 



and spin fluctuation scattering amplitude g^^. Here, n (w) is the Bose-Einstein distribution and Xi (Qi'^) is 
the dynamic susceptibility, given by Eq. 2 of the paper. Since the electric field defines a direction n, it is convenient 
to write the non-equilibrium distribution function as <i>k,A = ^k a'^'^' with /i = 1, 2, 3 denoting Cartesian coordinates. 

Given <&k it is straightforward to obtain the conductivity tensor cr^'' = ^ ( alTT ) ^k A*^k a- 

k,A 

To solve the Boltzmann equation (jSlI) . we write it in terms of differential operators (see, for instance, J. M. 
Ziman, Electrons and phonons: the theory of transport phenomena in solids) = (^Cimp + Csf^ ^'^j where X^^ = 

"^'^k A (^/k \/9£k.xj and $^ are vectors (/i = 1, 2, 3), but also matrices in the kA space. Ci is defined by the matrix 
elements: 



(«*'■),, = f E (*£,> 



K'.y 



k',A' 



k',A' 



-fO ^ j-A,A 

■'k,Aj '^i,k, 



(S3) 



In this language, it is convenient to introduce the scalar product ( X 



= Ek,A-^k,A^k.A: from which 



<i>''y. We consider that impurities are the dominant scatterers, while scattering by spin fluctuations 
perturbatively suppresses the conductivity. This is always true at low enough temperatures, where (T/7) is small 
compared to g^^^^- To lowest order, the solution of the Boltzmann equation is = -|-(5<I>'^, where <I>°''^ is the so- 
lution in the presence of impurities only and (5$^ is the solution of Ci^^S^^ = —C^i^'^'^. Therefore, the conductivity 
becomes: 
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imp 



(S4) 



The first term is the constant impurity contribution to the conductivity, Uinip, whereas the second term is the 
temperature dependent contribution from spin fluctuations (see Eq. 1 of the main text). The calculation of 
which satisfy — C'imp^^''^, is tedious but straightforward. We use the variational principle, as outlined in [2^, 
and the three-band model discussed in the paper, and obtain, for the electric field along x: 



0,x 

(2,3), 



V \ Stop + 2ff?,„p 



cos Q 



l + (5oT^2±252 cos 261 



(5?.„p + 2.g?^p) V Vl + ^o±^2cos 



26* 



COS0 



(S5) 



Here, B is the polar angle and VF^ £f, and v are, respectively, the Fermi velocity, Fermi energy, and density of states 
of the hole pocket. Si — Si/sp depends only on the band structure parameters and is given by: 



£F 



2fl 



1 - 



rrix + m 



2mxmy 



- ray 
^fhxfhy 



(S6) 



with mi = rrii/m and /i = [ijep- The residual conductivity due to impurities is: 



2 

imp 



i + 4(i + 5o) -2^(1 + ^0) -^i 



(S7) 



For an electric field parallel to the y direction, we can use the transformation properties of the non-equilibrium 
distribution functions under a 7r/2 rotation: = = and = with 9^9- 7r/2. 

We can now evaluate the matrix elements ( <I>"'^ Csf ) using the solutions in Eqs. IS5I After defining: 



(i) 



+ {2 + So ±62 cos 261') - 2 cos {9 - 9') y/l + So ± S2Cos29' 
±r] [cos 26* - 2 cos {9 + 9') 1 + So ±52Coa29' + cos 261' (l ± ^2 cos 261') 



(S8) 



we obtain (hereafter we omit the superscript in (f>j'g for convenience): 



/iMM _ r,. 

p pimp 
Pimp 



1=2,3 



with Pimp = o^im^p! t ~ (^-POo)^ and the pre-factor: 



(S9) 



5imp + 25i1 



l + 4(l + 5o) -2^(l + 5o) 



4^ 



(SIO) 



In Fig. 3 of the paper, we show the numerical solution of Eq. (IS9[) for the parameters discussed in the main text, 
using ^ [appearing in Eq. (S8)] that comes from the solution of the self-consistent equations (3) of the main text. 
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At low temperatures t <^1, we can obtain an analytical solution for the resistivity anisotropy. We consider that the 
system is in the orthorhombic paramagnetic phase, where the correlation length associated with each ordering vector 
is given by = T 'Pi where </? > is the Ising-nematic order parameter. To simplify our analysis, we consider a 
magnetic three-dimensional QCP. However, our results are more general and hold even away from the QCP. 

For the hard mode (short correlation length), ^^'^ = + ip, the t <^ 1 expansion of Eq. (|S9[) yields: 



y Pimp 



^imp 



hard 



C / dede' 



^3.6' 



f (3) , rjl 
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3/2 



(Sll) 



where we defined the auxiliary function Jq^\, — Lo^g\i 



(Ci/oo)" 



[see Eq. 

— 2 



(jSSP j. Thus, the contribution of the hard 
spin fluctuations to the resistivity vanishes as . For the soft mode ^2^^ = — lq the situation is different, since 
^2"^ -> as t ^ 0. Then, one has to divide the Fermi surface into two regions one around the hot spots Qhs^ 

f2) (2) 

with width t {fg g, ^ t), and another with fg g, 3> t, called the "cold region". Repeating the same steps that led to Eq. 
(|Slip . we find that the contribution from the cold region also vanishes as t^. Near the hot spots 9hs, an expansion 
similar to the one outlined in yields: 
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(S12) 



which is the Eq. (4) of the paper. Here, c > is a constant coming from the Jacobian of appropriate coordinate 
transformations. Since the contribution from the soft mode to the resistivity vanishes as 1?^"^ , while the other contri- 
butions vanish as t^, the former dominates the transport at low temperatures. Therefore, the resistivity anisotropy 
at t ^ 1 is: 
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with K > 0. Using Eqs. (jS5[) and the fact that the hot spots are at 9^^ = arccos((5o/'52)/2, we obtain: 



pVV _ pxx 
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(S14) 



Thus, for a fixed value of the mass anisotropy {S2), the sign of the resistivity anisotropy changes from positive to 
negative if the chemical potential {60) is negative and large enough. The "phase diagram" in Fig. 4 of the paper is 
obtained from this expression. 



